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Abstract 

Two different methods of obtaining "effective 2x2 Hamiltoni- 
ans" which include relativistic corrections to nonrelativistic calcula- 
tions are discussed. The standard Foldy-Wouthuysen transformation 
generates Hamiltonians which order by order in 1/M decouple the 
upper from the lower components. The upper left-hand block then 
defines an effective 2x2 Foldy-Wouthuysen Hamiltonian. In the sec- 
ond method the matrix element of the interaction Hamiltonian of the 
Dirac representation is evaluated between free positive-energy states 
and reduced to two-component form. The resulting expression (pos- 
sibly expanded in 1/M) then defines what we call the "direct Pauli 
reduction" effective 2x2 Hamiltonian. We wish to investigate under 
which circumstances the two approaches yield the same result. Using 
a generic interaction with harmonic time dependence we show that 
differences in the corresponding effective S-matrices do arise beyond 
first-order perturbation theory. We attribute them to the fact that 
the use of the direct reduction effective Hamiltonian involves the ad- 
ditional approximation of neglecting contributions from the negative- 
energy intermediate states, an approximation which is unnecessary 
in the Foldy-Wouthuysen case as there the 4x4 Hamiltonian does 
not connect positive- and negative-energy states. We conclude that 
at least in the cases where the relativistic Hamiltonian is known, us- 
ing the direct Pauli reduction effective Hamiltonian introduces spuri- 
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ous relativistic effects and therefore the Foldy-Wouthuysen reduction 
should be preferred. 

1 Introduction 

There are many cases in low- and intermediate-energy nuclear physics where 
one is confronted with the problem of incorporating relativistic effects in a 
nonrelativistic calculation. Such a situation arises, for example, when a rel- 
ativistic Hamiltonian description exists for one part of the interaction but 
another part can be realistically described only by a nonrelativistic Hamilto- 
nian. Naturally, as the energies involved in these processes increase, inclusion 
of relativistic corrections is expected to become more and more important. 

There are at least two conceivable ways to construct "effective 2x2 
Hamiltonians" in order to identify relativistic corrections to nonrelativistic 
Hamiltonians. 

If the full relativistic Hamiltonian, H, is known, the Foldy-Wouthuysen 
(FW) method [[[], ||, H f|, H, || |7|, |8|] provides a systematic procedure to block- 
diagonalize the Hamiltonian order by order in 1/M and hence to decouple 
positive- and negative-energy states to any desired order in 1/M. This is 
achieved by an application of successive unitary (in general time-dependent) 
transformations on the operator H — id/dt. This method has frequently been 
used in the context of electromagnetic processes in hadronic systems. 

Alternatively, one can think of evaluating the matrix element of the rela- 
tivistic interaction Hamiltonian between free positive-energy solutions. After 
reducing this matrix element to two-component form one may expand the 
resulting expression to any order in 1/M and use it as an effective Hamilto- 
nian Hj P to be evaluated between nonrelativistic Pauli wave functions. 
We shall refer to that procedure as the "direct Pauli reduction". 

The purpose of this paper is to address the question whether these two 
procedures yield identical results, and if not, which of the two should be 
preferred, as they both are used in the literature. For example, in electron 
scattering from nuclei, one evaluates matrix elements of electromagnetic op- 
erators ( "semi-de-relativized" according to the direct Pauli reduction proce- 
dure) between nonrelativistic solutions of the Schrodinger equation with a 
given NN potential ||, [H| [ll], [12], [131 , but there also exist treatments that 
make use of the Foldy-Wouthuysen procedure to obtain semi-relativistic op- 
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erators P, O, Similarly, in proton-proton bremsstrahlung the strong 



interaction part is conventionally treated using a nonrelativistic Lippmann- 
Schwinger approach, whereas for the electromagnetic interaction a semi- 
relativistic ansatz, involving either the Foldy-Wouthuysen method or the 
direct Pauli reduction method is made 0, [17], [L8|, [19|, |20fl . 

There are basically two observations that provide the key to our investiga- 
tion. The first observation is that there should be a clear distinction between 
Hamiltonians that are obtained starting from a relativistic one by means of 
unitary transformations (and which are still 4x4 operators) and what we shall 
refer to as effective 2x2 Hamiltonians (which are 2x2 operators, constructed 
as the upper left-hand parts of the corresponding 4x4 Hamiltonians). As 
we shall show, the direct Pauli reduction can be also viewed as an effective 
2x2 Hamiltonian obtained from a 4 x 4 unitary time-independent trans- 
formation of the relativistic Hamiltonian, by eliminating the contribution of 
the negative-energy sector of the relativistic Hamiltonian. This distinction 
extends to the S-matrices also, where the 4x4 Hamiltonians lead to the 
full relativistic S-matrix, whereas the effective 2x2 Hamiltonians lead to 
effective S-matrices. Such effective S-matrices, which are calculated in non- 
relativistic scattering theory, may, depending on the interaction, show effects 
of the neglect of the negative-energy sector. 

The second observation is that in most practical applications both these 
methods for obtaining effective Hamiltonians are not applied to the total 
(that is, electroweak-plus-strong) Hamiltonian but only to the electroweak 
part, in combination with a nonrelativistic approach to the strong inter- 
action. Physical observables, such as S-matrix, are then calculated in the 
framework of old-fashioned time-ordered perturbation theory. The examples 
presented above, namely, electron scattering from nuclei and proton-proton 
bremsstrahlung are clearly of that type. 

It is commonly believed that both types of effective Hamiltonians lead to 
identical effective S-matrices (and hence to identical observables) and thus 
that it makes no difference which is used in a nonrelativistic calculation. A 
justification for this assumption has been given in refs. [f2T], |22|. We discuss 
this assumption in quite some detail and find that in general it is not correct. 
In particular, in the cases where the full relativistic Hamiltonian is known, 
we show that the transformed 4x4 Hamiltonians all lead to the same S- 
matrix, which is (to the given order in 1/M) identical to the relativistic 
S-matrix. However, the effective S-matrices corresponding to the Foldy- 
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Wouthuysen and the direct Pauli reduction, which are those appropriate 
for a nonrelativistic calculation, differ from each other beyond first order 
in perturbation theory and only the Foldy-Wouthuysen effective S-matrix 
reproduces the full relativistic S-matrix to that particular order in 1/M. We 
discuss these issues in detail for the particular case of Compton scattering 



by a proton in ref. p3[ . In the cases where the full relativistic Hamiltonian 
is not known we cannot argue in favour of any of the different effective S- 
matrices in particular. In order to get a feeling about the importance of 
these differences, we compare the results of a proton-proton bremsstrahlung 
calculation using the two different methods. 

Our paper is organized as follows. In the next section we discuss gen- 
eral time-dependent unitary transformations. In this context we are mainly 
concerned with the Foldy-Wouthuysen transformation. In the third section 
we compare the first-order matrix elements of the effective Hamiltonians ob- 
tained with the Foldy-Wouthuysen and direct Pauli reduction procedures 
and investigate under which circumstances they lead to the same result. In 
the fourth section we argue about the need to go beyond first-order per- 
turbation theory. Using a generic Hamiltonian with the sole assumption of 
harmonic time dependence we explain how differences arise in the different 
effective S-matrices. In section 5 we give a general nonperturbative proof of 
why the S-matrix elements corresponding to the 4x4 transformed Hamil- 
tonians are identical with the full relativistic S-matrix. We conclude and 
summarize our main results in section 6. 



2 The Foldy-Wouthuysen Transformation 

In this section we discuss time-dependent unitary transformations with spe- 
cial emphasis on the Foldy-Wouthuysen transformation. 
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2.1 Unitary Transformations of the Schrodinger Equa- 
tion 

If we start with an equation of motion of the Schrodinger typeQ, 



Of 



H(t)Mt) >, (1) 



where we allow for an explicit time dependence of the Hamiltonian operator 
H(t), a unitary transformation T(t) 

|tt'(t)> = T(t)|*(t) >, 
T{t)T\t) = T\t)T(t) = 1, (2) 

will result in the new Schrodinger equation 

= ( r (W) Tt w - tT ^^r) >= > ■ 

(3) 

In many applications, such as e. g. the Foldy-Wouthuysen transformation 
[0, 0, [|, it is useful to parameterize T(t) as 

T(t) = e tS{t \ S(t) = S\t), (4) 

and then to expand the new Hamiltonian H'(t) in terms of S(t), 

H'(t) = H + t[S,H]-hs,[S,H}}- t -[S,[S,[S,H]}} + ... 

2 o 

+ £]] + ..., (5) 

2 o 

with 5 = dS/dt. 

As the unitary transformation T(t) may in general be time-dependent, 
one clearly finds 0, 121 



< *(t)\H(t)m) >^< ntwmnt) >, (6) 

which simply expresses the fact that if H (t) is the operator corresponding to a 
physical observable in the first representation, H'(t) is not the corresponding 

1 We regard the Dirac equation as a specific example of a Schrodinger type equation. 
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operator describing the same observable in the second representation [El], . 



In other words, as was pointed out by Nieto [pl|, it is the unitary transform of 
Hit), i. e. T(t)H(t)T^(t), which is physically equivalent to H(t), as it is this 
operator which yields the same matrix element between transformed states 
as H(t) between the original states. 

2.2 The Foldy-Wouthuysen Method 

The Foldy-Wouthuysen transformation |], ^ ||] provides a systematic method 
of finding a representation of the Dirac Hamiltonian in which positive- and 
negative-energy states are separately represented by two-component wave 
functions instead of the four-component wave functions in the ordinary Dirac 
representation. For the free Dirac equation^ 

l ^T~ = (C? ' P " + W* ^) = ff o*o(ar), (7) 

where at and (3 are the usual Dirac matrices, p is the momentum operator, 
and a; is a shorthand notation for (x, t), the transformation is exactly known 
and given by 

9™(x) = T * (x), 




To = \ ~^~ W E \ M . (8) 



where we have defined the operator E = y/p 2 + M 2 . In the new Foldy- 
Wouthuysen representation the free Dirac equation has the simple form 

= PyJp + w^V) = h% w < w {x). (9) 

As (3 is of block-diagonal form, eq. (|S|) is just the direct sum of two Hamil- 
tonians ±\/p 2 + M 2 . The positive/negative-energy solutions of eq. @ are 
of the form 

*r <+, w = ( x<+ o w ) ! 



2 States and Hamiltonians without specific superscript labels will always be assumed to 
be in the Dirac representation. 
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where X ( x ) an d are two-component spinors. Note that Hq W is simply 
given by TqHqTq because T is time-independent. 

If the Dirac Hamiltonian contains an explicitly time-dependent inter- 
action Hj(t), the transformation T will depend on time, T = T(t), and, in 
general, a closed form for the transformation T(t) leading to a block-diagonal 
form of the Hamiltonian is not known. However, Foldy and Wouthuysen [|TJ 
developed a systematic procedure to construct a new Hamiltonian which is 
block-diagonal to any desired order in 1/M. The idea is to split the Dirac 
Hamiltonian into its odd and even components^, 

H(t) =H + H^t) =(3M + Oit) + £j(t), (11) 

where 0(t) = + Gi(t) = a ■ p + Oi(t). It is then assumed, that the 
interaction potentials Oi and Si do not contain powers of 1/M that are 
smaller than (1/M)°. Furthermore, the interaction has to be weak enough 
in the sense that the magnitude of each of the time and space Fourier com- 
ponents of the interaction potential is considerably smaller than the mass 
of the nucleon Jl|. It is then understood that we mean by [1/M] "terms of 
order E re f/M" , where E re f is some reference energy smaller than M. The 
procedure consists of first applying the transformation 

r«(t) = e isW ^\ 

SV)(t) = -iP^4, (12) 

to eq. (|IID, the result of which is then written as 

h fw{i) (t) =(3M + qW ( t) + g (i) ( t ) _ ( 13 ) 

Using eq. (^) it can be easily shown [III |^, that the transformation of 
eq. (0) is constructed such that the odd operator of eq. (p~3[ ) is of order 1/M. 
The procedure is then repeated with a new transformation T^ 2 \t) which is 
exactly of the same form as eq. flT2p except that in ^( 2 )(t) the new odd 
operator 0^(t) of eq. © appears. As each successive application reduces 
the leading power of the odd operator of the resulting Hamiltonian by one 

3 Odd operators O, such as e. g. cti, couple large and small components whereas even 
operators £ (e. g. (3) do not. The following identities are useful in the derivation of the 
Foldy- Wouthuysen transformation: [O, 0] = 2(9/3, [£, 0\ = 0. 
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unit, after applying this method n times one obtains a Hamiltonian which is 
block-diagonal to order 1/M n_1 . For example, after four transformations the 
Hamiltonian which does not contain any odd operators up to and including 
order 1/M 3 reads 

= ,( M+ ^_|l) +£ ,__L [0 , [0 , £/] 

where Oi stands for dOj/dt. 

3 First— Order Matrix Elements 

It is in general not possible to explicitly solve the Dirac equation with an 
arbitrary time-dependent potential. On the other hand in many cases it is 
sufficient to construct a nonrelativistic approximation including lowest-order 
relativistic corrections. Such an approximation is often formulated in terms 
of a 2 x 2 nonrelativistic effective Hamiltonian H e ^ which is to be used in 
nonrelativistic perturbation theory. One way of constructing such an effec- 
tive Hamiltonian to be used in low- and intermediate-energy applications is 
to evaluate the matrix element of the potential operator between positive- 
energy solutions of the free Dirac equation, perform a two-component reduc- 
tion and then interpret the result as the matrix element of an effective 2x2 
Hamiltonian between nonrelativistic (positive-energy) states (direct Pauli 
reduction). Commonly the so constructed effective Hamiltonian is expanded 
in a series in 1/M and only the terms to a particular order are kept. 

3.1 Reduction of the First— Order Matrix Element 

To be specific we consider the expression 

J<Px*l f (2,t)H I {t)* 0i &t) = J d 3 x^p^ t ) ToHl (t)T^F w (x,t), (15) 

where ^f 0i and are positive-energy solutions of the free Dirac equation 
and \?o$ and are the corresponding solutions of the Foldy-Wouthuysen 
Hamiltonian (see eq. (||) and (§)). Expressions of the type of eq. (|I5]) typically 
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appear as building blocks of a perturbative treatment of the S-matrix. From 



eq. (15) we find that instead of explicitly reducing the matrix element, we 



may as well look at the operator Hf = T Q Hj(t)TQ. This operator, when taken 
between the states anc ^ ^o/^j leads to an effective 2x2 Hamiltonian 

fftff-P which j s j us t w hat one would get by performing a two-component 
Pauli reduction on the left-hand side of eq. (|15D. 

For a general interaction, Hj(t) = Oj(t) + Sj(t), a reduction through 
order 1/M 3 yields 

Hf(t) = Sj + ^iOM-^p^O^S^-^iOliOM}} 
' -pi [Oq, Oi\\ + [1/M 4 ] + odd terms. (16) 



16M 3 ' 

Here the terms written explicitly are all even operators. The odd terms 
do not contribute when eq. ( |16"D is evaluated between free positive-energy 
states of the Foldy-Wouthuysen representation. However, as it is only these 
states which are considered in a nonrelativistic calculation, one can regard 
the upper left-hand block of Hf{t) as an effective Hamiltonian H e /^~ P {t) to 
be used with two-component Pauli spinors. 

3.2 Comparison with the Foldy— Wouthuysen Hamilto- 
nian 



The Foldy-Wouthuysen Hamiltonian of eq. fll4|) may be written as 

H FW(4) = H FW + H F W{t) + jjFW^ + H F W{t) + R FW ^ + [1/M ^ 

where the subscripts % indicate the order in which the interaction potentials 
Oj and Si appear in the corresponding Hamiltonian. The linear interaction 
Hamiltonian is given by 

H[ w (t) = S I + ^{O ,O I }-^ r2 [O ,[O ,S I ]+ i d I ] 

ft r/n2 f/n /-on , n / n/r^ 



;{Oi{O 0l Oj}}+[l/M% (18) 

OlVl 

lparing eq. ( JIBD and eq. i 
Hf{t) differ by 



8M 3 

Comparing eq. (TO) and eq. flTBD we see that the two operators H[ w (t) and 
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= -g^[Oo, Oj] + j^fsPl [Ob, Oi]] + odd terms + [1/M 4 ] 

= ([C [Oo, <?/]] - i[0o, 0/]) + odd terms + [1/M 4 ]. 

(19) 

In obtaining the last line of eq. ([03) we expanded eq. (0) and made use of 
0,1 P 2 - 

It has often (tacitly) been assumed that AH\(t) of eq. ( |T9| ) vanishes when 
evaluated between free positive-energy Foldy-Wouthuysen states. In which 
case, in a nonrelativistic calculation to first order in the interaction, it should 
make no difference which of the two Hamiltonians Hf(t) and H[ w (t) is used 
for the interaction. We will now investigate under which circumstances this 
assumption is correct. 

Most applications of the Foldy-Wouthuysen transformation are concerned 
with either static (atomic physics, hydrogen atom) or harmonically time- 
dependent external potentials. If we assume the potential 0r(t) to have a 
harmonic time dependence, we find for the time derivative Oj = ^iuOi, 
where the upper/lower sign refers to the absorption/emission of a quantum 
of energy to. In this case the matrix element of AHi(t) reads 

JptEf-EkTv) J d"x^Y\x 1 t)[0 Ql O I \^J{x,t) + [l/M% 

(20) 

From eq. (|20|) it is seen, that for a potential with harmonic time dependence, 
both Hamiltonians will yield the same first-order matrix elements if the 
energies of the initial and final state match with the absorption/emission 
of a quantum of energy u, or to put it in a somewhat sloppy way, if energy 
is "conserved" in the transition, or alternatively if the commutator [Oq, Oj] 
vanishes, which in general we do not expect to happen. This observation has 
also been made in ref. 0. 

In second-order (old-fashioned) perturbation theory, energy is not con- 
served at a vertex |26|. j27fl . This clearly suggests, that in practical second- 
order calculations, such as e. g. proton-proton bremsstrahlung, it may make 
a difference whether one uses the linear effective Foldy-Wouthuysen Hamil- 
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tonian Hi^~ FW (t) deduced from eq. ( |T8| ) or the Hamiltonian H e /^~ P {t) ob- 
tained from eq. (|l6l). 

Our result seems to be somewhat different from the earlier work of Nieto 



[f21|j and of Goldman [fj2j, who appear to have shown that the results from 
the two first-order Hamiltonians are the same, i. e. that the matrix element 
of AHi(t) = in general. To understand this difference, consider again this 
matrix element. Using the last part of eq. (|19"D one can write 

J d^x^^AH^t^^t) = 

' d'x ^P(x, t) ([A, H FW ] + <A) , t) + [1/M% (21) 

where to the order being considered A = [Oq, Oi]/8M 2 . A more general 



form for A is given in ref. |[22|| , but is unnecessary for the purposes of this 
discussion. 



The authors of both ref. PI| and ref. make the assumption that 

[A, H F w ] + *A = 0, (22) 

which they refer to as the "Schrodinger statement" , and thus conclude that 
< ^J f w/t (t)|A J ff 1 (t)|^^ iy (t) >= and consequently that it makes no differ- 
ence which of the Hamiltonians is used. 

It is clear however from eq. (19) that eq. ( p2j) cannot be correct in general 



as an operator equation. If one considers it instead as a matrix equation, 
then the Schrodinger equation can be used to replace the Hamiltonian by a 
time derivative and one obtains 

<*™(t)\([A,H FW ]+i*)\< W (t) >= 
. d 
'It 



% T + < *$F{t)\A\*gr(t) > +[1/M 4 ]. (23) 



Again, however, the right-hand side of this equation is not identically zero 
in general. If one assumes a harmonic time dependence of the interaction, 
as we have done above, then the time dependence of the matrix element 
is just e l< - £ ' / ~' EiTW - ,i and the time derivative simply brings down the energy 
factor i(Ef — Ei =f to). This leads exactly to our previous eq. (^) and to 
our previous conclusion that the matrix element of AH\ vanishes, and the 
results using different Hamiltonians are equivalent, only when the transition 
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conserves energy. In the more general case when the states are off-energy- 
shell, as will be the case when such matrix elements are used in second-order 
perturbation theory, different Hamiltonians may be expected to give different 
results. 

As a specific example of a case where the choice of Hamiltonian does 
make a difference, let us consider the interaction of a proton with an external 
electromagnetic field, 

Bj[t) = e$(x,t)--^(3a-B(x,t), 

Oj(t) = -ea-A(x,t)+i^3a-E(x,t), (24) 

where e > is the elementary charge, k = 1.79 the anomalous magnetic 

moment in units of a nuclear magneton, and E = — V$ — A and B = V x A. 
Using this specific example the even part of AHi(t) is given to order 1/M 3 
by 

eB - -> -> 

AH x {t) = -j^[p 2 ,-iV ■ A + ia-px A-ia- Axp\ 

B{p- E + E -p + a -V x E) 



16M 3 

ie , 

_ 8M^ 



-iV ■ A + io -p x A- id ■ A x p). (25) 



In all modern proton-proton bremsstrahlung calculations 0, [IT], 0, [19|, ^0 
the electromagnetic part of the interaction is obtained by reducing the rel- 
ativistic nucleon-nucleon-gamma vertex to obtain an effective 2x2 Hamil- 
tonian which is used with two-component spinors in a nonrelativistic cal- 
culation. In some cases, e. g. ref. |l6 l, the effective Hamiltonian ff^f~ FW 



obtained via a Foldy-Wouthuysen reduction as in eq. (|18| ) was used and in 
others, e. g. the H e /^ P obtained via a direct two-component reduction 



as in eq. ( |TED was used. These two should differ by the AHi of eq. ([05]) or 
eq. dH). 

In figure 1 the results of such a bremsstrahlung calculation are shown 
for a particular kinematic situation, corresponding to a coplanar equal-angle 
geometry with outgoing protons measured at 03 = O4 = 6° and with an in- 
cident laboratory energy T LAB = 280 MeV. The calculations of ref. [TH] have 



been used taking a) the Foldy-Wouthuysen Hamiltonian for the electromag- 
netic interaction of the proton or b) the effective Hamiltonian of eq. (|i~6|). As 
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can be seen from eq. (p0| ) the difference between the two calculations is of 
order 1 / M 2 multiplied by the amount by which energy is not conserved at the 
electromagnetic vertex. For the kinematics chosen in figure 1 the difference 
is of order 7 %. 

We conclude that as long as one considers only first-order matrix ele- 
ments, such as those appearing in fi + p — > n + v or electron scattering 
from an on-shell nucleon [24] it does not make a difference which reduc- 
tion scheme is used. This is true because in a perturbative treatment of the 
S-matrix the time integration enforces energy conservation at the vertex of 
lowest-order perturbation theory. However, as we have seen in figure 1 the 
situation is different for higher-order processes such as e. g. proton-proton 
bremsstrahlung. 

One might properly fault this particular numerical example because the 
strong interaction is not treated consistently with the electromagnetic inter- 
action, i. e. not as a reduction of some relativistic interaction. Thus it is 
important to look at more detail at second-order processes in which both 
interactions can be treated in the same way. That we do in general in the 
following section and for the specific case of Compton scattering by a proton 
in ref . fE3[ . 



4 Second— Order Perturbation Theory 

In the previous sections we described how to obtain various effective interac- 
tion Hamiltonians for use in nonrelativistic calculations and showed that the 
first-order matrix elements of these Hamiltonians differed only in situations 
when the states involved did not conserve energy. Such a situation arises in 
second- (or higher-) order old-fashioned time-ordered perturbation theory 
where such first-order matrix elements connect to an intermediate state. Of 
course, the total energy is conserved for each (higher-order) diagram, but en- 



ergy is not conserved at the individual vertices of such a diagram |26|, \Ff\ . It 
is this fact which may lead to different predictions for the effective S-matrix 
and thus it is important to see what happens in such second-order processes. 

There are two questions which have to be addressed. Firstly, which kind 
of transformation of the 4x4 Hamiltonian leads to the same relativistic 
S-matrix? As long as we are dealing with Hamiltonians originating from a 
unitary transformation of the states (see section 2.1) one would expect that 
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the full relativistic S-matrix would be invariant under such transformations. 
That is indeed the case as we will show in second-order perturbation theory 
in the first part of this section and by more formal manipulations valid to all 
orders in sec. |^. 

This invariance of the full relativistic S-matrix has been observed by oth- 
ers, and is consistent with our usual understanding of quantum mechanics. 
However it has led to the extremely misleading, and incorrect, assumption 
that the physical observables calculated with various effective Hamiltonians 
will be the same and thus that it makes no difference which effective Hamil- 
tonian is used. 

Hence we are led to the second question which is of great practical im- 
portance. Do different effective 2x2 Hamiltonians lead to the same non- 
relativistic effective S-matrix? We show in the last part of this section that 
in fact this is not the case. The reason is that the construction of effective 
Hamiltonians is not always done solely via a unitary transformation of the 
states. Thus it does make a difference which effective Hamiltonian is used in 
nonrelativistic calculations. Only the Foldy-Wouthuysen Hamiltonian (and 
other block-diagonal Hamiltonians obtained from it) reproduces the full rel- 
ativistic S-matrix to a given order in 1/M. Thus it is only such Hamiltonians 
which can be presumed to give correct physical results in nonrelativistic cal- 
culations. 

4.1 Full Relativistic S— Matrix in Second— Order Per- 
turbation Theory 

We now restrict our consideration to second-order perturbation theory. Thus 
we assume that the relativistic 4x4 Dirac Hamiltonian is given as before 
by H = H + Hj(t), where H is the free Dirac Hamiltonian. To be definite 
we will focus on processes involving an incoming quantum of energy u a and 
an outgoing one of energy Ub- Thus Hi(t) will have a term ~ H a e~ luJat and 
one ~ Hbe lu)bt a.s well as others as required for Hermiticity. Formally we could 
construct Hj{t) in second quantized form and use appropriate creation and 
destruction operators to pick out the pieces needed. Instead we will treat 
Hj(t) somewhat loosely with the understanding that, at the end, by second- 
order we will mean that we will discard all terms except those containing one 
power of the coupling associated with H a and one power of that associated 
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with Hi,. It may be that H a = Hb, as for example in Compton scattering 
where H a = Hi, ~ H^nn- Alternatively H a may be different from H^ as in 
pion photoproduction where we could take H a ~ H lNN and H b ~ H wNN . 

With these assumptions the full relativistic S-matrix is given through 
second order in the interaction in the standard fashion, e. g. ref. @, eq. (6.57), 
as 

S fi = Sfi-iftfyVlfMHiMVoiiy) 

-i J d^z^^HMSriy - s)/^*) (26) 

where \l/o/ and \l/oi are positiv e-energy eigenstates of the free Dirac equation 
and where Sf is the free Feynman propagator. To proceed further we need 
to examine how the states, Hamiltonian, and propagator transform under 
the unitary transformation T(t). 

We define the states ^/q p (x) to be eigenstates of the free Dirac Hamiltonian 
Hq. They may have positive or negative energy and can be written as 

Vfffrt) = e^$g } (f), (27) 

where E p = +y / p' 2 + m' 2 and 

(2k) 



where ■u/ + - ) = u and ' = v with v)u = v^v = 1 are the usual Dirac 
spinorsf]. Spins will always be summed, so all explicit spin dependence will 
be suppressed. These states satisfy a completeness relation 

spins 

Under the unitary transformation T , which is all that will be needed, we 
have ^>Q p W<<± \x) = T ^f^\x). This changes only the spinors which become, 
using the specific T of eq. (^j), 

u(p)-fXj and «(p)-(°). (3°) 



Note that we use a different normalization convention in comparison with ref. @. 
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with x a two-component spinor and a two-component null spinor. We will 
always assume that To is obtained from T(t) by turning off the interaction, 
which thus removes the time dependence. For most of this section however it 
is not necessary to specify T or T(t) explicitly. Thus we will use \P' and H' 
for example for states and Hamiltonians obtained via general transformations 
Ol. We will reserve ^> FW and H FW for situations where it is 



as m sec. 



important that the states have the form eq. (|30D and that the Hamiltonian is 
block-diagonal to some order in 1/M, and thus where we require the specific 



transformation T of eq. (§) and the specific T(t) of sec. 

To determine how the propagator S F transforms it is easiest to first 
expand in the complete set of states just defined. Thus starting with, 
e. g. eq. (6.48) of ref. 0, and using a standard expansion for 6(t — t') we 
obtain 

■■'> ni -n^i,-i ) V I c -iEJt v ~t z ) ®Qp\y)®Qp^ ( z ) 

Po + it 

(31) 



S F {y-z)P = —J dp^pe'^-^ 



spins 

-po + ie 



We now define 

S' F (V - z)P = T S F (y - z)f3Tl (32) 

It is clear by applying T . . . Tq to the right-hand side of eq. (|3~ID that 
S' F (x — y)(3 has exactly the same form as S F {x — y)/3, only with the states 
$ replaced by $' and so is the correct free Feynman propagator in the $' 
representation. 

Now consider the Hamiltonian. To allow an expansion to second order 
we write the transformation variously as 

T(t) = e iS = T e iS = r e^+ & +-), (33) 

where S n is Hermitian and of order n in the interaction. Starting from 
H' = T{t)H{t)T\t) - iT{t)t\t) and inverting, we find for H 

H = e~ iS { T^H% + ^ lS -^~ \ e iS - (34) 
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We now write H' as H' = Hq + H[ + H 2 + . . . , where H' n is of order n in the 
interaction and expand if in a fashion analogous to eq. @ to obtain through 
second order 

E 1 = {t^H[T -AH[} + {t^H' 2 T ~AH 2 }, (35) 
AH[ = T^AH^ =i[3 h Ho] - & , (36) 

AH' 2 = i [S 2 ,H }-S 2 - i -[S 1 ,AH[]+t{S 1 ,T^H[T ]. (37) 

H' 2 is second-order in the interaction and is a contact term which is generated 
by the transformation. Hi is however first-order in the interaction, so the 
second bracketed expression in eq. (B5f) must be zero. Hence AH' 2 = T ] Q H' 2 T . 

The next step is to substitute this result for Hi into eq. (|26|) and use 
eq. ( |3"^) for S' F (x — y)f3 to obtain a result for that part of the full relativistic 
S-matrix which is second-order in the interaction, 

% = -if d^^ f {y)H'M^{y) 

-i J d A yd A z^ f {y)H' x {y)S' F {y - z^H'^z^z) 

+AS fh (38) 

where 

A% = i J d 4 y^l f (y)AH 2 (y)^ 0l (y) 

-i J d 4 yd±z¥ 0f (y){ - T^H[(y)T S F (y - z){3AH[(z) 
-AH[(y)S F (y-z)pr$H[(y)T 

+AH[{y)S F {y - z)AH[{z)}^ 0i {z), (39) 

Observe now that if ASfi = 0, as we will in fact show, then the right-hand 
side of eq. ( |3"8"D is exactly the S-matrix that one would compute to second- 
order perturbation theory using the transformed Hamiltonian H' and the 
transformed states. Comparison with eq. (^) shows then that the S-matrix 
is invariant under the unitary transformations of the type we have considered. 

Thus consider ASfi. It is sufficient to look only at the AH x AH term 
as the others are all similar. This term becomes after using the explicit 
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expression for Sf, eq. (|31|), and extracting the time dependence from the 
interaction and the states 



/ dt y dt z dpod 3 px 

E 



—iA+t v +iB+t : 



spins 



Po + ie 
i(-) ^/ nsHi 



— e M 

-Po + it 



+ c. t. 



-2ni5{E f + u b - Ei 



E 



d 3 p x 
(+)i 



spins 



+ 



E f + u b - E p + ie 
Ef + uj b + E p - ie 



+ c. t. 



(40) 



where v4± = p — Ef — uj b ±E p and £>± = po — Ei—u a ±Ep. Here c. t. indicates 
that there is a standard cross term which has not been written explicitly. 

Using the time dependence of Si, which is linear in Hi, the matrix ele- 
ments of AH[ can be evaluated as 



<$$\AH[ a \<S> 0i > 



< ®$\i[Si a ,H } +iu a S la \$oi > 

i{Ei + u a TE p ) < ®$\S la \$ i >, (41) 



< $ of \AH[ b \$ { f > = < $ of \i[S lb , H ] - iuj b S lb \$ ( f > 

= - l (E f + uj b TEp)<% f \S lb \^ f > . (42) 

These matrix elements are analogous to that of eq. ( p0[) but, since they 
connect to intermediate states $op\ energy is not conserved and the energy 
factors do not vanish. 

Finally, using these equations in a symmetric way (so as to end up with 
the commutator) and using the completeness relation eq. (P9f), one finds for 
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eq. 



- 2iri5(E f + Ub -E i -u a )<$ of \-^[S 1 , AH[]\<5> 0l > . (43) 

To simplify the notation, the formula has been rewritten with the full inter- 
action, thus reintroducing the H a x H a and H b x H b terms which should in 
the end be dropped. 

The other terms are evaluated in exactly the same way so that one obtains 

ASfi = 2m6(E f + u b -E i -u a )<$ of \AH' 2 + l[S 1 ,AHtt 

-t[S u T^H[T ]\^ 0l > 
= 2m5(E f + u b - Et- uj a ) < $ Q f\i[S 2 , H ] - i(u b - uj a )S 2 \&oi > 
= 2n5(E f + u b - Ei - w a )(E f + uj b -E l - u a ) < % f \S 2 \%i > 
= 0. (44) 

Here the first reduction uses the definition of AH' 2 of eq. flT?|) and the second 
uses overall conservation of energy. 

In summary, we have in this section considered the full relativistic S- 
matrix in second-order perturbation theory (though to all orders in 1/M). 
For definiteness we focussed only on those processes involving an incoming 
quantum of energy u a and outgoing quantum of energy uj b . In that context 
we could show explicitly that the AS/i of eq. (J38|) was zero. The remaining 
terms on the right-hand side of that equation give just the full relativistic 
S-matrix expressed in terms of the transformed Hamiltonian H\ and the 
transformed states to second order. Thus we have shown explicitly that 
the full relativistic S-matrix is unchanged under these transformations in 
second-order perturbation theory, which is a special case of the more general, 
but more formal, proof to be given below. 

Note that nowhere in this proof did we use any of the properties of 
T(t), except those of eq. ([33|). Thus the result that the full relativistic S- 
matrix is invariant holds for the transformation T(t) of sec. |2.2| which block- 



diagonalizes Hi and produces the Foldy-Wouthuysen Hamiltonian Hf w . It 
also holds for T(t) = T , where the specific T of eq. (H) is meant, which 
produces Hf . 

One other observation is worth making. The matrix element of AH is 
essentially the difference between the matrix elements of Hi in the $o basis 



19 



and Hj in the <3>' basis. From the details of the proof that ASfi = 0, 
as outlined in eqs. ( [TOD - d44|) , one can see that these AH terms cancel 
the energy denominators and thus generate contact terms. In general these 
contact terms which make up AH 2 came both from positive- and negative- 
energy intermediate states. This is somewhat in contrast to the generally 
held belief that it is the negative-energy intermediate states alone which are 
responsible for contact terms. It agrees however with the specific result we 



have obtained in ref. [p3|1 . There the specific example of Compton scattering 
by a proton is described in detail and we can also see at which order in 1/M 
the various terms arise. 



4.2 Comparison of Effective 

Nonrelativistic S— matrices in Second— Order Per- 
turbation Theory 

We now want to return to the main question, namely whether or not different 
reductions of the relativistic interaction to effective nonrelativistic interac- 
tions lead in practical calculations to different results. To do this we need 
first to review how these effective Hamiltonians are obtained and used. 

Usually one starts with a supposedly known relativistic interaction and 
makes a transformation on it followed by a projection onto the upper left- 
hand block which leads to a 2 x 2 effective interaction Hamiltonian, Hj , 
which is then used together with positive-energy two-component states to 
construct an effective S-matrix, Sf{, according to the usual rules of non- 
relativistic scattering theory. Terms of higher than leading order in 1/M 
included in the effective Hamiltonian provide relativistic corrections to the 
lowest order nonrelativistic result and one hopes that if enough terms are in- 
cluded S e JJ will approach the full relativistic, and thus presumedly correct, 
result, Sfi. 

Thus the question which is relevant for practical calculations is whether 
or not these effective S-matrices are the same for different reductions of the 
relativistic Hamiltonian, i. e. , for different Hj, and whether they equal the 
full relativistic S/j. 

In the Foldy-Wouthuysen case the new Hf w is obtained via a trans- 
formation T(t) of the type we have been considering, constrained by the 
condition that to some order in 1/M the 4x4 Hamiltonian Hf w is block- 
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diagonal. The 2x2 effective Hamiltonian }{ e /^~ FW is then just the upper 
left-hand block, i. e. that part which has nonzero matrix element between 
free positive-energy states § FW . Note that by virtue of the fact that T(t) 
depends on the interaction, ff^f~ FW w ill have terms of higher order in the 
interaction, even though Hj was first-order. 

An alternative approach often used is to make a direct Pauli reduction 
of the relativistic matrix element < $ /|#i|$oi >=< ® F f W l^o#/7o >. 
Thus in this case the 4 x 4 Hamiltonian Hf = TqH^ and the 2 x 2 Hj ff ~ P 
is just the upper left-hand block of this. Note that in this case Hf is not 
block-diagonal and that H e /^~ P may or may not be expanded (or truncated) 
in powers of 1/M. Furthermore H e /^~ P is linear in the interaction, so there 
are no contact terms generated naturally. 

To see how the effective S-matrices calculated from these effective Hamil- 
tonians are related to the full relativistic S-matrix we start with eq. (B3) 
above with ASfi = and use eq. (^l|) and eq. (]32| ) to expand the propagator 



into positive- and negative-energy states. The time dependence can then be 
extracted and the formula simplified in exactly the same way the AH x AH 
term was treated. The result for the second-order contributions only is 

S ft = -2mS(E f + u b -E i -u a )S.<^ F f w \H' 2 \^ pw > 
htJ V E t + uj h - E n + ie 



spins 



E f + u b - E p + ie 

v > 

+ c. t. 



Ef + u b + E p - ie 

" (45) 

The term with H 2 and the terms with & F p W ^ involve only positive-energy 
states, which have only upper components. Hence these terms pick out just 
the upper left-hand block of H[ and H 2 and consequently make up just the 
effective S-matrix one would calculate using the 2 x 2 H e / f . Thus we can 
write 



S fi = S e P -2m8{E f + uj b -E i - 



OJa) X 



spins 



Ef + u b + E p - ie 
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(46) 



Now it is obvious that for the Foldy-Wouthuysen reduction the last term 
vanishes, since there H[ w is block-diagonal and since a block-diagonal Hf w 
does not connect positive- and negative-energy states. Thus in this case 
S fi = Sf- FW to some order in 1/M. Hence the Foldy-Wouthuysen proce- 
dure gives the correct effective S-matrix when }j e /^~ FW is used cLS 9i 2 X 2 
nonrelativistic Hamiltonian in standard nonrelativistic scattering theory. 

Note that the only thing needed to get this result is an interaction Hf w 
which is block-diagonal to a given order in 1/M. Thus a further transfor- 
mation on the Hamiltonian which preserves the block-diagonal nature to the 
same order in 1/M will also give the same correct S-matrix S e ^J~ FW . Thus 
it is clear that the Foldy-Wouthuysen algorithm is not unique. It has been 
observed by Barnhill and others (see ||] and references therein) that the or- 
der by which different, noncommuting odd operators are eliminated from the 
relativistic Hamiltonian may lead to such a "freedom" in the choice of the 
FW Hamiltonian. 

In contrast to the Foldy-Wouthuysen reduction, the direct Pauli reduc- 
tion leads to a 4 x 4 Hamiltonian Hf = T HjTq which is not block-diagonal 
in general. Thus the terms in Sfi involving negative-energy states will not 
vanish. Thus if we use Hj consistently in a nonrelativistic scattering 
theory we will get an effective S-matrix Sf{ which is not equal to the 
(correct) relativistic result and which is therefore different from that ob- 
tained using Hjff~ FW , by virtue of the ^ terms in eq. (jiff ). In other 
words, using H e /^~ P in a nonrelativistic approach amounts to omitting the 
negative-frequency contribution of the covariant calculation. One may ex- 
pect that this approximation will lead to unacceptable results in situations 
where the original interaction Hamiltonian in the Dirac representation pro- 
duces a significant coupling between upper and lower components. 

We conclude from this that in cases in which both interactions are known 
relativistically, so that one can do a Foldy-Wouthuysen reduction on the 
complete interaction, as for example Compton scattering, use of Hj 
gives the correct result to given order in 1/M, the same result as would be 
obtained in the complete relativistic theory. In contrast the Pauli reduction 
effective Hamiltonian Hj*~ P gives different, and incorrect, results. 

In many practical cases however one knows only one of the interactions 
relativistically. Thus for example in proton-proton bremsstrahlung or in 
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the electrodisintegration of nuclei the electromagnetic interaction is known 
relativistically but the strong interaction is not. In such cases the reduction 
is done on the known interaction alone, and one cannot generate the contact 
terms which involve both interactions. Thus it appears that for these cases 
one cannot prove rigorously that one effective Hamiltonian is more correct 
than another 5 . 

To see how this works in a still somewhat simplified case, suppose that 
the interaction part of the Dirac Hamiltonian is given as before by harmonic 
interactions of the form Hi = H\ a + Hn,. Then by carrying out exactly 
the same steps used to evaluate the AH x AH term in eq. fl4*0| ) we find for 
the part of the full relativistic S-matrix which is first-order in each of the 
interactions 



S 



fi 



-2ni5(E f + u b - Ei 



d 3 p x 



E 

spins 



< <hf\Hn\<$>$ >< $$\H la \$oi > 



E f + uj b - E p + ie 



i(-) 



< $ /|#lfc|$0p ; >< Spp^-HlalSoi > 



+ 



Ef + ujh + E p 



ie 



+ c. t. 



(47) 



We suppose that we do not know the interaction Hu in a relativis- 
tic form. This means that we must approximate the full matrix elements 
< c&o/l-^ibl^op' 1 > an d < &of\Hib\$op^ > in some way. Typically one would 
neglect the latter and evaluate the former by using a nonrelativistic matrix 
element calculated, say, from a potential. Then one would reduce the known 
interaction if la to an effective Hamiltonian. 

Using eqs. Q33) and (|55|) one has H ia = TqH^Tq in the Pauli reduction 
case and H la = V H™T - AH[ a with AH[ a = i[S la , H ]-S la in the Foldy- 
Wouthuysen case. Substituting these into eq. ( fT7D above and using eq. ((4T|) 
to evaluate the AH terms we find for the part in brackets 



< $ 



- Oi 



> 



< $ 



0/ 



Ef + uo h 



E v + ie 



5 Of course there may be other guiding principles, such as gauge invariance, which put 
a constraint on the form of the effective Hamiltonian. 
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(48) 



The result using; Ht f J~ p is of the same form, except that the second term 
does not appear and H^J~ FW is replaced by H.\^J~ P . 

The first term in this expression is what one would usually use in second- 
order perturbation theory in cases where only one of the interactions is 
known. It would lead to an approximate S-matrix, s^ Pprox , which differs 
from the full relativistic S-matrix by contact terms. Using an expansion 
analogous to eq. (H) and explicit forms for the leading parts of S\ a and T 
it is straightforward to show that the leading term in T iSi a TQ is odd and 
of order Oi/M whereas the first even term is down by one power of 1/M. 
Thus the contact term involving the negative-energy intermediate states is 
likely to be the most important, whereas the one involving positive-energy 
intermediate states (which appears only for H^J~ FW ) is down by 1/M. One 
can see this explicitly in Compton scattering [p3 |. 

Above we showed that when both interactions are known fails to 

give the correct effective S-matrix because of the neglect of the terms involv- 
ing negative-energy intermediate states. On the other hand H e /^~ FW gives 
the correct answer, since it is constructed as the upper left-hand block of a 
Hamiltonian which does not connect positive- and negative-energy states, 
and since it contains the contact terms explicitly. Now when only one of the 
interactions is known the natural way of calculating S^ Pprox in second-order 
perturbation theory forces one to neglect these contributions for both effec- 
tive Hamiltonians, and leads to an error which to leading order is the same 
in the two cases. 

In most practical calculations, particularly when the unknown interaction 
is the strong interaction, second-order perturbation theory is not sufficient. 
In such cases, an example being proton-proton bremsstrahlung, one must 
sum the strong interaction to all orders. Thus one normally would replace 
< $ j|if lfe |$Qp > in eq. ( f4~8|) by a full T-matrix. The cross term gives the 
interaction in the other order and there is an additional term, called double 
scattering in proton-proton bremsstrahlung, which has a T-matrix element 
on both sides of the matrix element of H\J . An alternative, but essentially 
equivalent approach which is often used is simply to take matrix elements of 
H\U between eigenstates of the free-plus-strong Hamiltonian. 

Neither of these approaches will in general give the correct (relativistic) 
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answer, as long as the effective operator H^J is obtained as a reduction 
of the relativistic Hamiltonian between free states, because of inconsistency 
of the operator and the wave functions. In other words, approximations 
will be involved in getting from the full relativistic S-matrix to the effective 
nonrelativistic one. Regardless of the choice of H^l, contact terms will be 
missed. 

It appears then that here, in contrast to the situation when both interac- 
tions are known relativistically, we cannot argue rigorously that }j e /^~ FW is 
any more correct than H e /^~ P . At best it may be preferred simply because it 
is no worse an approximation when the interaction is incompletely known and 
is in principle the correct choice when the interaction is completely known. 



5 Transformation of the S— Matrix under a 
Unitary Transformation 

In the previous section we have explicitly shown to second-order perturbation 
theory that the S-matrix remains invariant under a unitary transformation 
of the states^. Here we will show this in general without making use of a 
perturbative approach. We start with the Dirac equation 

d^(x) 

i—^ = (H + Hj(x, 0)*(z) = H(x, 0*(x), (49) 

where Hq = a ■ p + /3M, x = (x°,x). We control the switching on 
and off of the interaction in the remote past and distant future through 
Hi(x,^) = Hj(x)e~^,^ > 0. The defining equation for the full propagator 
in the presence of interaction reads |2], |3| 

- Hj(x, o) S FtHl (x, y)(3 = 5\x - y), (50) 




6 In fact we can even allow for more general transformations (see also ref. (£9|). It 
turns out that the following proof is still correct if we demand T(t) to be an invertible 
transformation with the property of eq. (|57|). Note however that To remains a time- 
independent unitary transformation. One then only has to replace T* by T" 1 in the 
proof. 
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and a perturbative solution satisfying the Feynman-Stiickelberg boundary 
condition is given by 

S F>Hl {x,y) = S F (x - y) + J d A x 1 S F (x - x 1 )(3H I (x 1 ,£)S F (x 1 - y) 

+ J d 4 xid A X2S F (x — X2)(3Hi(x2,£,)S F (x2 — xi)[3Hj(x\, £)S F (xi — y) 
+ ..., (51) 

where S F (x — y) is the free Feynman propagator [Q, |3|. For < x° < y® the 
formal solution for may be constructed as 

*(ar)=i J d 3 y 1 S F>Hl (x,y 1 )(3^+\y 1 )-i J d 3 y 2 S F!Hl (x,y 2 )(3¥-\y 2 ), (52) 

where ^f^(yi) and ^^(y 2 ) are superpositions of positive- and negative- 
energy solutions of the free Dirac equation, respectively If we specify the 
boundary conditions as 

lim ¥ + \x) = ty ( t\x), 

x°—>—oo 

lim ¥-\x) = 0, (53) 

cc°— >oo 

where ( x ) * s a positive-energy eigenfunction of the free Dirac equation 
with eigenvalue Ei > 0, we can either describe the scattering of a particle or 
pair annihilation in the potential, 

Sfi = lim lim / d 3 x ^q^\x)^(x) for particle scattering, 

£— >0+ x°— >oo J 

Sfi = lim lim / d 3 x tyi* (x)ifj(x) for pair annihilation, (54) 

with ^qj^ a negative-energy solution of the free Dirac equation. Likewise, 
specifying the boundary conditions 

lim ¥~\x) = ^ { ~\x), 

x u — >oo 

lim ¥ + \x) = 0, (55) 
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leads to a description of the scattering of an antiparticle or pair creation in 
the potential, 

S n = lim lim /d'x*H*(x)*(x) for antiparticle scattering, 
Sfi = lim lim / d 3 x ^f r ^\x)'9(x) for pair creation. (56) 

When defining the S-matrix elements of eq. (|54]) and (|56|) , it is assumed that 
the system was in an eigenstate of H Q with positive energy in the remote 
past or will be in an eigenstate of H with negative energy in the distant 
future, respectively. It interacts with the potential at intermediate times and 
and evolves under the influence of Hq in the distant future or remote past, 
respectively. The S-matrix is the mathematical idealization of extending the 
interaction over the complete t-axis [28|. However, it is important to realize 
that the limit £ — > + has to be taken at the end (see e. g. ref. 0, p 165 f). 

If we introduce a unitary transformation T(x) which depends on the in- 
teraction Hamiltonian Hi(x,£), or parts of it, with the property 

lim T(x) = T , (57) 

where T is a time-independent unitary transformation, we will show that 

S'fi = S fi , (58) 

provided the boundary conditions and S^ are defined in an analogous way 
as in eq. (|53|) - (p6|), e. g. for particle scattering 



lim ^ ,(+) (x) = *& +) (o0, 

x°^—oo 

lim #4-) (3) = 0, 

S'fi = lim lim [ d 3 x^p\x)V'{x). (59) 

1 i— >0+ x°— >oo J 1 

In eq. ( |59"D ^'^(x) and ( x ) are positive-energy eigenstates of H = 



Of 

TqHqTq and ty'(x) satisfies the Dirac equation 

= (H' + H' T (x, 0W(x) = H'(x, 0*'(x), (60) 
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with H'(x,£) = T(x)H(x,£)T'(x) — iT(x)dT>(x)/dt. The extension to the 
other cases is straightforward. We then find for the S-matrix element 5^(0 

= lim [ d 3 x^\x)T^T(xMx). (61) 

J x°— >oo J J 

For any arbitrarily small but finite £ eq. (|57|) leads to 

lim T f T(x) = 1 (62) 



x°— >oo 



and we therefore obtain 

^(0 = Sfiit). (63) 

Taking the limit £ — > + in eq. (|63|) yields the desired result, eq. (|58|) . 

The above derivation is extremely simple, but it does not explicitly reveal 
the importance of H'(x, £) of eq. (|60|) , in particular, it does not show why it is 
H'(x, £) and not T(x)H(x,£)T\x) which enters the calculation of the primed 
S-matrix element. For that reason we provide a slightly more complicated 
derivation which, however, gives more insight into the role played by H'(x,£). 
We once again only discuss the case of particle scattering. Using eq. (|2|) and 
the boundary condition of eq. (|59]) we may write S 1 ^ as 



S' fi = • lim lim [d 3 xd 3 y^\x)S' m (x,y)^\x) 

= i lim lim [d 3 xd 3 y^\x)T^S' F ^(x,y)PT ¥ t\x), (64) 

where S' FH ^(x,y) is defined by an equation analogous to eq. (|5UD, namely 

(i^ -H>- Hfe, 0) S F , H ,(x, y)(3 = 5\x - y). (65) 

Clearly, at this point one realizes that H'(x) rather than T(x)H(x,£)T*(x) 
is the relevant operator in the defining equation for the full propagator. It 
may easily be shown that 

S'vH'fav)/ 3 = T ( x ) S F,H I (x,y)/3T\y) (66) 
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solves the equation of motion for S' FH t , eq. (|65|), provided Sf,Hj satisfies 
eq. (ISO). We then insert eq. (EH) into eq. (El) to obtain 



,(+), 



Oi 



X). 



S' fi = zlim lim [ d 3 xd 3 y¥+ ) \x)TZT(x)S F!Hl (x,y)(3T\y)T y[ 

(67) 

Once again we make use of eq. fl57|) for any arbitrarily small but finite £, 

lim T*(y)T = 1 = lim rjr^), (68) 

and obtain 

S' f M) = * lim [ {Px<Py*$\x)S F , Hj (x,y)l3*£\x) 

y°— oo 

= %(0, (69) 

which is identical with eq. (|63|). 

In conclusion, we have provided two nonperturbative arguments for the 
S-matrices to be equal provided the states are related by a unitary transfor- 
mation. 



6 Summary and Conclusions 

We have addressed in this paper the question of how one obtains effective 
nonrelativistic Hamiltonians from known, or partially known, relativistic in- 
teractions. Such Hamiltonians are important because they allow the use of 
nonrelativistic formalisms, and make connection with interactions which may 
be known only in some nonrelativistic approximation, and yet at the same 
time incorporate relativistic corrections to some order in 1/M. 

We discussed two different ways of obtaining such effective 2x2 Hamil- 
tonians. The first is the Foldy-Wouthuysen interaction which is obtained by 
a unitary time-dependent transformation of the states. This approach de- 
couples the upper from the lower components order by order in 1/M and so, 
by projecting out the upper left-hand block, leads naturally to an effective 
Hamiltonian H e /^~ FW to be used with two-component positive-energy wave 
functions. 
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The second method consists of making a two-component Pauli reduction 
of the matrix element of the interaction Hamiltonian in the Dirac representa- 
tion between free positive-energy spinors. The result is then regarded as an 
effective Hamiltonian Hj^~ p to be used between nonrelativistic Pauli wave 
functions. This approach is equivalent to making a unitary transformation of 
the states using the time-independent transformation T and then projecting 
out the upper left-hand block to get H e /^~ F . Usually the result is expanded 
in 1/M, typically to order 1/M 2 or 1/M 3 , though in some applications this 
step is omitted, i. e. the complete expression without a 1/M expansion is 
used. Both of these methods for obtaining effective Hamiltonians and intro- 
ducing relativistic corrections have been used for a variety of processes in the 
literature. 

It has been tacitly assumed in the literature that both approaches yield 
the same results, i. e. that it makes no difference which of the two effective 
Hamiltonians is used. We have seen that this is in general not true and that 
in fact it does make a difference which is used, at least in the case when the 
initial relativistic interaction is known completely. 

For first-order matrix elements differences in the results obtained using 
the two different effective Hamiltonians arise only when at least one state 
is off-energy-shell, or in other words only when energy is not conserved, 
as the difference of the matrix elements calculated with different effective 
Hamiltonians is proportional to the sum of initial energies minus the sum of 
the final energies at a vertex. 

However, in higher-than-first-order, time-ordered perturbation theory, 
energy is not conserved at each individual vertex. Furthermore intermediate 
negative-energy states contribute and so we might expect effects from the 
projection eliminating the effects of the negative-energy states which is used 
in getting the effective Hamiltonians from the relativistic ones. We found 
in second-order perturbation theory that the effective S-matrices, obtained 
by using the 2x2 effective Hamiltonians in usual nonrelativistic scattering 
theory were in fact different. The result obtained using }j^f~ FW agreed with 
the full relativistic S-matrix while that obtained with Hj did not. 

At first glance this seems inconsistent with the perception that the ob- 
servables should not change under a unitary transformation. As we showed, 
both in second-order perturbation theory and via a general argument, the full 
relativistic S-matrix is invariant under such time-dependent unitary trans- 
formations, of which T(t) leading to Hf w and T leading to Hf are examples. 
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This apparent conflict was resolved via an understanding of the additional 
approximations necessary in going from the full relativistic S-matrix to the 
effective S-matrix calculated nonrelativistically. To get S e ^J one must ne- 
glect the negative-energy intermediate states. Since Hf w is block-diagonal 
and does not connect positive- to negative-energy states anyway, this ap- 
proximation has no effect and in this case the effective S-matrix w is 
the same as the full relativistic one. In contrast Hf does connect negative- 
and positive-energy states. Thus the additional approximation needed to get 
Hj^~ P , namely the projection for positive energies which picks out just the 
upper left-hand block, leads to the neglect of some nonzero terms and thus 
to an effective S-matrix Sf{ which is different from the full relativistic 
S-matrix. 

The results just described apply only to cases where the full relativistic 
interaction is known, as for example Compton scattering or pion photopro- 
duction. In many practical situations only part of the relativistic Hamiltonian 
is known. Typically, such a situation may involve a nonrelativistic potential 
approach for the strong interaction and at the same time a relativistic treat- 
ment of the electromagnetic interaction (including the anomalous magnetic 
moment and possibly on-shell form factors). 

In such cases the situation is ambiguous. Lack of knowledge of one of 
the interactions makes it impossible to calculate the contact terms explicitly 
and hence they must be neglected whichever effective Hamiltonian is used. 
Use of H e /f~ FW requires one to drop an additional contact term, which is 
however of order 1/M times the main contact terms which are dropped in 
all cases. For some processes gauge invariance allows one to construct non- 
unique contact terms which can be added in by hand. However in general 
there may be numerical differences when different effective Hamiltonians are 
used. Unfortunately it does not seem to be possible to determine in a rigorous 
way which is best for these situations, though there may be a philosophical 
predilection to use the Foldy-Wouthuysen approach even in the absence of 
full information as that is the correct approach when the interaction is fully 
known. 

In conclusion then, our analysis indicates that when the interaction is 
fully known relativistically use of the Foldy-Wouthuysen reduction in order 
to incorporate relativistic corrections to a nonrelativistic treatment leads 
to correct results. In contrast using the matrix element of the interaction 
Hamiltonian in the Dirac representation between free positive-energy states 
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as an effective nonrelativistic Hamiltonian gives incorrect results. However 
when only part of the interaction is known, as is the case in many practical 
situations, both methods miss the leading order contact terms. Thus in this 
situation one cannot argue rigorously that one should be preferred over the 
other. 
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Figure 1: Cross section and analyzing power for proton-proton 
bremsstrahlung for an incident laboratory energy of 280 MeV and a coplanar 
equal angle geometry for the outgoing protons, plotted versus the photon 
angle. The calculations were done using the full Bonn potential and the for- 
malism of ref. except that Coulomb corrections are not included. The 
solid (dotted) curve corresponds to using n e ^~ FW (Hj^~ P ) for the electro- 
magnetic interaction. Note the suppressed zeros on both axes. 
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